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I. INTRODUCTION 



We present a study of the Casimir effect^"- occuring when an imperfect Bose gas fills 
a region between two infinite parallel plane walls separated by distance D. In particular, 
we are interested in the properties of the fluctuation-induced force between the walls as 
function of distance D, and of the thermodynamic state of the imperfect Bose gas in the 
vicinity of the Bose-Einstein condensation point. Our work is thus a contribution to the 
theory of critical Casimir effects in condensed matter physics. 

The case of a perfect Bose gas has been analyzed by Martin and Zagrebnov^. The authors 
considered the grand canonical potential fi(T, L, D, ji) of a perfect Bose gas filling a box 
(rectangular parallelepiped) with linear dimensions L x L x D. The temperature and the 
chemical potential are denoted by T and /i, respectively. Periodic, Dirichlet, and Neumann 
boundary conditions were considered. The main result of their study was the determination 
of the large D asymptotics of the grand canonical potential per unit wall area 

U)(T,D,/j) = hm (1) 

L—>oo L 

It turns out that Bose-Einstein condensation induces a qualitative change in the behavior 
of the excess grand canonical free energy u s (T, D, /i) 

u a (T,D,ti = [u{T,D,y)-Du) b (T,y)\ , (2) 

where Ub denotes the potential Q per unit volume evaluated in the thermodynamic limit. 
The excess grand canonical free energy u s (T, D, /j,) contains contributions due to each single 
wall confining the system and the interaction between the confining walls. 
In the one-phase region corresponding to \i < 0, the force 

F(T,D,,) = - d ^ T gD D -^ (3) 

decays exponentially fast for D — > oo, whereas in the presence of condensate, when fi = 0, 

there appears a long-range Casimir force 

F(T,D,0) = _A 

k B T D 3 U 

The Casimir force as defined in Eq.Q is evaluated per unit area and has dimension of 
pressure. The amplitude A depends on the type of boundary conditions imposed at the 



walls and is otherwise universal. For example, in the case of periodic boundary conditions 
the amplitude A is negative and corresponds to attractive Casimir forces. The analysis 
presented in 8 was extended in 9 , where the dependence of the Casimir force on thickness D 
and temperature was presented with the help of universal scaling function of variable D/£, 
where £ is the bulk correlation length of density fluctuations. This extension contains both 
the Casimir effect at criticality and its asymptotic behavior off criticality as special cases. 

Our object in the present study is to generalize the above analysis to the case of an 
imperfect (mean-field) Bose gas. The repulsive pair interaction between identical bosons 
is described in this model by associating with each pair of particles the same mean energy 
(a/V), where a > 0, and V denotes the volume occupied by the system. The Hamiltonian 
of the imperfect Bose gas 15 composed of N particles 



H = H + H mf (5) 



is the sum of the kinetic energy 



k 

and the term representing the mean-field approximation to the interparticle interaction 

E mf = ^\ ■ (7) 

The symbols {n^} denote the particle number operators and the summation is over one- 
particle states {k}. We follow here the generally adopted definition of the imperfect Bose 
gas where the exact number of pairs N(N — l)/2 is replaced in H m f by N 2 /2 (see the 
corresponding comment in^). As far as the applicability of the mean-field approximation 
is concerned, an important rigorous result is to be noted. It turns out that the mean field 
theory can be equivalently formulated with the help of a repulsive binary potential subject 
to the so-called Kac's scaling V 7 (r) = 7 3 t> (jr), where v(r) is the actual binary potential. 
When the dimensionless parameter 7 approaches zero, the potential V 7 (r) becomes weak and 
long range, leading exactly to the mean field theory based on the Hamiltonian in Eq.([51) with 
constant a equal to the 7-independent space integral of the potential Vy— . This remarkable 
connection clarifies the physical content of the mean field model studied here. 

The thermodynamics of the imperfect Bose gas has been rigorously studied showing 
that Bose-Einstein condensation persists in the presence of mean- field interaction 13 . The 



corresponding critical density n c and critical temperature T c are the same as in the perfect 
Only the critical chemical potential is shifted from /x = in the perfect gas 
case to a positive value \i c = an c . We can thus expect the appearance of Casimir forces 
much as in the case of the perfect gas. However, the imperfect Bose gas is a qualitatively 
different system which motivates the present study. Indeed, the presence of the term H m f 
makes the system's Hamiltonian superstable which, in turn, implies equivalence between 
canonical and grand-canonical ensembles^ii. This is not the case for the perfect gas where 
the grand-canonical description leads to unphysical macroscopic fluctuations in the overall 
density in the ground state, whereas the canonical ensemble predicts physically satisfactory 
behavior—*^. Moreover, the thermodynamics of the imperfect Bose gas is defined for all 
values of the chemical potential, whereas the states of an ideal Bose gas are confined to the 
region /j < 0. In the present paper we thus investigate Casimir forces in physically more 
satisfactory case of interacting mean-field Bose gas. 

In Section II we present an original derivation of the bulk grand-canonical free energy 
density Ub(T,fi) using the steepest descent method. The same method is then applied to 
determine the Casimir force and study its scaling properties in Section III. The last Section 
IV contains concluding comments. 



II. BULK PROPERTIES OF THE IMPERFECT BOSE GAS VIA THE 
STEEPEST DESCENT METHOD 

We consider the imperfect Bose gas composed of N identical particles of mass m contained 
in L x L x D rectangular parallelepiped of volume V = L 2 D; D is the distance between two 
L x L walls. 

Owing to the structure of the Hamiltonian (jSJ) the grand canonical partition function 
H(T, L, D,fi) can be conveniently written in the occupation numbers representation 

oo 

~(T,L,D,(j,) = Vexp(/3/iiV)y; exp 

where ^/„ k \ denotes summation over all sets {^k} of occupation numbers of one-particle 
states {k} under the constraint N = ^ k ^k- The kinetic energy in state k equals e k = 
h 2 \k\ 2 /2m = h 2 k 2 /2m, and & = l/k B T. 
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-/3 ^2 n ^ ek + 



2 V 



(8) 



l(T,L,D,fi) = exp ( y-/^ 2 ) exp 



The partition function can be rewritten as 

T Rn V ii _1 

Z (T,L,D,N) , (9) 
where Z (T, L, D, N) denotes the canonical partition function of a pefect Bose gas 

(10) 



Z (T,L,D,N) = J2 {nk} exp[-P^2 



It is useful to rewrite Eq.(j9]) by introducing an arbitrary real parameter a in the following 

way 



i(T,L,D,fi) = exp 



2a 



(/i — a)' 



(11) 



x 



zZ ex P \ ~ 



N=0 



(3a 
2V 



V 

N (u-oc) 

a 



exp(PNa)Z {T,L,D,N) 



where the rhs of Eq. ffTTj) actually does not depend on a. Such a representation of 
H(T, L, D, fx) has been already used in the proof of condensation taking place in the imperfect 
Bose gas 13 . 
The identity 

„2\ Ifo /-x 



exp — 



r 

48 



dq exp(—5q + i^q) 



(12) 



valid for arbitrary 5 > 0, permits one to establish relation between H(T, L, D, fi) and the 
grand partition function of a perfect gas S (T, L, D, fi) 



S (T, L, D,fi)=J2 exp(/3/iiV) Z (T, L, D, N) 



(13) 



N=0 



evaluated for a complex value of the chemical potential (a + iq). Indeed, the following 
sequence of equalities holds 



(T, L, D, fi) = exp 

V/3 



[3V 



00 \\r a roo 

x, V 2: '" ■ 



2a 
dq exp 



(fi-a) 



—q 2 + iq/3 [N-—(fi- a)} + (3aN 
2a a 



(14) 



Z (T,L,D,N) 



exp 



W 2 

M 

2a P 



27ra 



c?g exp 



— (iq + a) -(iq + a) 

2a a 



E Q (T,L,D,a + iq) 
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-%) exp 



2^ 



XL 

27ca 



o+ioo 



dt exp 



E {T,L,D,t) 



Relation (JHj) between H(T, L,D,fi) and H (T, L,D,t) will be the basis for subsequent 
analysis. The parameter a appears here only in the contour of integration. Due to its 
arbitrariness we can choose a < to ensure the convergence of the series defining the 
analytic continuation of the perfect gas partition function ([131) . 

We will apply here the method of steepest descent (saddle point method) to determine 
the large L and large D asymptotics of the partition function S(T, L,D,/i). This, in turn, 
will yield the bulk free energy density (per unit volume) as well as the excess grand canonical 
free energy density (per unit area). 

To illustrate the method we begin by evaluating in this section the bulk free energy 
density. This is conveniently done by considering a cubic volume V = L 3 . We choose the 
coordinate system with axes oriented along the edges of the cube adopting in each direction 
periodic boundary conditions. The rhs of Eq.( fl4l) can be then rewritten in the following way 



2a ^ 



V 



2vr/3a 



/ ds exp[-V ip b (s)} , (15) 

J 8a—ioo 



S(T, L,L,/x) = — i exp 
where 

1 / s 2 \ 1 1 



Here A = h/ y/2iTmkBT denotes the thermal de Broglie wavelength , and g K (z) is the Bose 
function 



oo u 
Z K 



k=l 

Equation ffTB 7 ]) follows directly from the formula 

1 11 

— In S (T,L,L,fi) = —g 5/2 [exp(0n)] - — ln[l - exp(^i)] (18) 

valid for V/X 3 ^> 1 (see e.g.— >^). It is well known that the theory of Bose-Einstein conden- 
sation in a perfect gas requires separate examination of the contribution [— ln[l — exp(s)]/V] 
stemming from the one-particle ground state. 



In the steepest descent method^ one has to find the saddle point so at which the derivative 
of the analytic function (pb(s) vanishes 

o = viM = JL (- So + M - l i93/2M ] - i . (19) 

This equation can be rewritten with the help the critical density n c for the perfect Bose gas 

A 3 n c = £3/2(1) = C(3/2) = 2.612... (20) 

in the following way 

so + J_ = ^3/2[exp(g )] + _\_ j" exp(so) j ^ 



exp(s J 



exp(s J 



f3an c an c £3/2(1) V n c 

Depending on the value of parameter (fi/an c ) one has to consider separately two cases: 
(fj,/an c ) < 1, and (///an c ) > 1. 

In the first case, when (/i/an c ) < 1, there exists a unique real negative solution (s < 0) 
satisfying the equation 

so + = ff3/2[exp(s )] ^ 



/3an c an c £3/2(1 
This follows directly from the analysis of the point of intersection of the straight line 

f(s ) = + (23) 

pan c a n c 



f v £3/ 2 [exp(g )] 

j(so) = "^MlT (24) 



with the curve 



represented in Fig.l. Clearly, when so < 0, the second term on the rhs of (12 ip does not 
contribute in the bulk limit V — > 00. We also note that the saddle point s (T,fi) tends to 
zero from below when approaching the imperfect gas condensation point fi = an c . 

Applying the method of steepest descent^ to determine the large V asymptotics of the 
integral in Eq. (115j) we find that the bulk free energy density 

uj b (T, fi) = - lim — k B T In E(T, L, L, //) = -p (25) 



is given by 

2a (3 2 

where p denotes the pressure and the superscript < indicates that we consider the case of 
[L < an c . 



a" 1 f s 2 \ 1 

ufiT, fji) = -— - — I -± - SoPnj - — £ 5/2 [exp(s )] , (26) 




FIG. 1. Graphical solution of Eq. (|2ip . Two straight lines represent the lhs of Eq. (|22p in two 
cases: the dotted line corresponds to = 0.9 and f3an c = 0.5 (one-phase region); the dashed line 
corresponds to = 1.1 and f3an c = 0.5 (two-phase region). The continuous line corresponds to 
the rhs of Eq. (]2ip in the limit V — > oo. The density of the condensate fi/a — n c is proportional to 
the segment denoted by the bracket on vertical axes. 



The bulk density us^ (T, jj) can be equivalently rewritten as 

Ub{T,p) = -p K {T,y) = -^x 2 (T,/i) - j^g 5/2 [exp(P(fj,-ax(T,fj))] , (27) 

where x(T, /x) satisfies the equation 

x = ^3 s- 3 /2[exp[/3(/i - ax)} . (28) 

In fact, a straightforward calculation shows that 

|^<(T, M ) = -*(T, M ) (29) 

so that x(T, n) = n(T, fi) is the number density of the gas expressed as function of chemical 
potential and temperature. 
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In the complementary case (n/an c ) > 1 the solution of Eg. (12 II) is obtained by taking the 
limit so — > 0, L — >• oo , and requiring the equality 



- = n c + 
a 



lim -Jr- 

sq— ¥0,L— >oo Ll 



exp(s J 



n c + lim — 

so— >0,L— >oo L 



- exp(s )_ 
n c + (density of condensate) 



(30) 



1 - exp(so) 

see the dashed line in Fig.l. The precise meaning of the double limit So — > 0, L — > oo 
defining the density of condensate in the perfect gas theory can be found 

The bulk free energy density in the region (fi/an c ) > 1 is thus obtained by taking the 
limit so — > in Eq.()26l) , so that 

■ 2 1 



or 



(3X- 



■C(5/2) 



(31) 



(32) 



where Po(T, /x = 0) is the perfect gas pressure in the presence of condensate. Eq. (I3~2~j) implies 
the relation n — fi/a, valid in the presence of condensate. 

The above results for the bulk density of the grand canonical potential yield equations 
of state (I27|) . (l28l) . (l32|) . in full agreement with rigorous results derived by more complex 
methods (see e.g.—). 



III. CASIMIR FORCE BETWEEN PARALLEL PLANE WALLS 

In order to evaluate the Casimir force (per unit area) acting between infinite parallel 
plane walls separated by distance D one has to determine the corresponding excess grand 
canonical free energy density (per unit area) defined by 

~Q(T, L, D, /i)' 



o; s (T, D 1 \x) = lim 

L— >oo 



Du b {T,fi) . (33) 



V- 

The system surface properties, and the excess grand canonical free energy in particular, 
depend on the boundary conditions imposed on the confining walls. Here we concentrate 
on the case of periodic boundary conditions although the generalization to the Dirichlet 
and Neumann cases is straightforward. We choose the coordinate system whose z-axis is 
perpendicular to the confining walls. 



It is convenient to rewrite the partition function (|T4|) in the following way 

Q(T,L,D,fiY 



E(T,L,D,n) = exp 



knT 



(34) 



-% exp 



(3L 2 D 
2a 



l 2 d r? a+ ' 

2TT/3a 



ds exp(— L 2 D (p(s)) 



where 



(p(s)= <p b (s) 



- - 

A 3 I D 



92 t exp ( s ~ P e ( k z))) - ^5/2[exp(s)] 



(35) 



+ 7TH E < ln t 1 - ex P ( s " ^))] - ln [l " ex P( s )]> 



L 2 D 

Here e(fc 2 ) = h 2 k 2 /2m, and fc z = 2nn/D, n = 0, ±1,±2, •■• due to periodic boundary 
conditions. Equation ( 135]) follows directly from the equality 

\nE (T,L,D,s/f3) = - £ ln{l - exp[s - (3e(k z )}} - (j) <7 2 [exp[s - /3e(^)r36) 

fc z ^o k z ^ ' 

valid for L/A ^> 1. The first term on the rhs of (1361) corresponds to the k x = k y = 

contribution to the series defining the logarithm of the perfect gas partition function 



In ~ (T, L, D, s/P) = -J2 M 1 ~ ex P(* " 



(37) 



Finally, the bulk density (fb{s) which we added and subtracted when writing ( 1351) has been 
defined in Eq. (Tl6l) . 

Taking in Eq. fl35l) the L — > oo limit we obtain the formula 



<p(s) = p b {s) 



- - 

A 3 I D 



92 [exp[s - 0e(k z ))) - # 5/2 [exp(s)] 



(38) 



corresponding to infinite walls. With the help of the Jacobi identity^ 1 ^ 

oo ^ oo 

exp(— nun 2 ) = —= exp(— nn 2 / k) 

n=— oo 

valid for k > 0, and relation 



A 

— I -nn" 
D 



Eq.flHB]) can be rewritten in the following form 



oo oo 1 



A 3 



exp I sr 



r=l n=l 



D 2 nn 
X 2 r 



(39) 



(40) 



(41) 
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The saddle point s of function <p(s) satisfies the equation 



2 ~ ~ i / D 2 7rn 2 \ 

T— 1 « — 1 ' >■ * 



(42) 

r=l n=l x ' 

In fact, we are interested in the behavior of s for large values of the dimensionless width 
D/X, and the sole property of s we will need is the fact that in the limit D/X — > oo the 
point s approaches the bulk saddle point s . Indeed, within the steepest descent method the 
solution s is inserted into (p(s) in order to evaluate the surface contribution fl33|) to the grand 
canonical potential. In the case of periodic boundary conditions the excess grand canonical 
free energy density oj s (T, D, p) represents exclusively the fluctuation-induced interaction be- 
tween the walls; there is no contribution to u s (T, D, p) stemming from two noninteracting 
wall-boson gas interfaces. 

The above remarks combined with Eq. (|4ip imply that for large D/X the dominant contri- 
bution to cu s (T, D, fx) has the form 



u} s (T,/i,D) 



2 D k B T 
A 1 



oo oo 1 



exp 



D 2 nn 2 



s r - 



X 2 r 



(43) 



r=l n=l 

where So is the bulk saddle point defined by Eq. (ll9j) . From the previous analysis we know 
that for n/an c < 1, so takes a negative value ensuring the convergence of the series in 

For D/X ^> 1 the sum over r on the rhs of Eq. (|43|) can be transformed as^ 1 ^ 



oo ^ 

Y] — exp 

r=l r ~ 2 



D 2 7rn 2 



s r - 



X 2 r 



dp— 

i) ps 



exp 



2n- 



D 



'per 



exp ( — 2n- 



SoP - 
D 



D 2 Tin 
~}?p~ 



21 



(44) 



'per 



where 



A 



K 



per 



(45) 



a/ vt|s | 

is the characteristic decay length. 

We have already noticed that sq — > when the condensation point p = an c is approached 
from below. It is instructive to determine the structure of the the solution of Eq.(l22l) in the 
vicinity of the condensation point. With the help of the asymptotic relation— 



23/2 (exp(so)) « -2^ {-so) 1 ' 2 +C(3/2) 



(46) 
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we find from Eq. (f22|) that 



1/2 _ (an c -fi) C(3/2) 

s — o 1 /2 ' l^'J 



Near the condensation point the behavior of the decay length ( 1451) is thus described by 

_ an c 27T 1 / 2 

K per = A • 48 

(cm c - /i) C(3/2) 

Thus upon approaching the condensation point from the one-phase region the decay length 
K per diverges like (an c — /J,) . We notice that the value of the relevant critical exponent 
vimp = 1 differs from that corresponding to the perfect Bose gas, where n per ~ (— /i) -1 / 2 
with the corresponding critical exponent vp = 1/2 (see 8 -). 
For the imperfect Bose gas the Casimir force 

D 



du<(T,D^) d II ^ 1 + 2nD/ Kper 
^ = dD = ksT dD 1 ^ 22 e ^ 

n=l 



-2n- 



K-per _ 



(49) 



can be rewritten with the help of the scaling function 9 



= -J2 1 + 2 3 nX exp(-2mr) , (50) 



7m° 

n=l 



peri 



such that ^(0) = -C(3)/vr, ^'(0) = 0, and V(x > 1) » -2j-Kxe- 2x . Putting ac = D/« 
we find in the one-phase region (normal phase) near the condensation point 

F{T k ° T ^ = ±[29(x)-x9'(*)) ■ (51) 

We note that the scaling function 2^(x) — x^'(x) in Eq. (l51j) depends on the ratio x = 
D I ' Kper of the width D and the decay length n per , see Eq.f l45|) . An interesting problem is to 
determine the correlation length £ characterizing the decay of the density-density correlation 
function of the imperfect Bose gas and then to express the variable x as function of D/£. 
This would allow one to calculate the universal scaling function for the Casimir forced. 

On the other hand, in the two-phase region (in the presence of condensate) one observes 
the power-law decay 

F(T,D^>an c ) _2((3) 1 

which is exactly the same, and with the same amplitude (universal when divided by kpT) 
as in the perfect Bose gas 8 at /i = 0. For comparisons with analytical and numerical results 
for the Casimir amplitude and the scaling function for the XY universality class which 
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have been obtained via various methods, inter alia via Monte Carlo simulations and the 
field-theoretical renormalization-group theory for O(N) symmetric systems, see^ 1 ^— . 

The above analysis can be straightforwardly extended to the case of Dirichlet and Neu- 
mann boundary conditions. We thus present here only the final results. In these cases - 
contrary to the case of periodic boundary conditions - one obtains nonzero values for the 
single wall-Bose gas surface free energy density cr(T, /x). The excess grand canonical free en- 
ergy density is the sum of 2<r(T, /x), and the contribution due to effective interaction between 
the walls. 

For a < (T, /x) one finds 

^(T,^) = ±^-£ 2 [exp(s )] , (53) 

where sq is the solution of Eq. (122j) . and the (+) and (— ) signs correspond to the Dirichlet 
and Neumann boundary conditions, respectively. The expressions ( 1531) for the surface free 
energy density should be compared with the corresponding results obtained for the perfect 
Bose gas, see Chapter 5 in^. One concludes that although in general these expressions are 
different, they become identical when evaluated in the two-phase region, i.e., fi = for the 
perfect gas, and /x > an c for the imperfect gas. 

The contribution to the excess grand canonical free energy density due to the effective 
interaction between the walls is given by 

k B T 1 + 2uD/k d n 
-^D^L ~3 ~ exp(-2rx J D/K AJV ) , (54) 

n=l 

where the expression for the decay length kd.n is the same in the Dirichlet and Neumann 

cases 

K>D,N = 7=== , (55) 

2v 7r l s o| 

and differes by factor 2 from the corresponding expression K per in the periodic case (see 
Eq.gHD). 

In the presence of condensed phase, where /x > /x c = an c , one observes both in the 
Dirichlet and in the Neumann case the same power law decay of the corresponding Casimir 
force 

F(T,D,ti>an c ) C(3) 1 
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with the same amplitude (universal when divided by kpT) as in the perfect Bose gas for 
these boundary conditions (see 8 ). 

IV. CONCLUDING COMMENTS 

The aim of our study was to investigate Casimir forces induced by fluctuations of an 
imperfect Bose gas contained between two infinite parallel plane walls situated at distance 
D, and compare the results with those obtained previously in Ref.- for a perfect gas. We 
based our approach on exact relation (11411 between the grand canonical partition function 
of an imperfect (mean-field) gas, and that of an perfect gas. The derivation of Eq. (1141) was 
based on the fact that the mean-field Hamiltonian, Eq.Q, allowed the use of the occupation 
number representation (for a systematic description of the so-called diagonal models see the 
review^) . 

We first applied to Eq. (jT4l the steepest descent method to determine the bulk free energy 
density. In this way we obtained the mean-field equation of state, previously established 
rigorously by other methods. It might be interesting to note that the self-consistent relation 
( |28l) has been derived and studied for an interacting Bose gas in a much more general 
context 25 . 

The evaluation of the excess grand canonical free energy density under periodic boundary 
conditions, presented in Section III, used in an essential way the fact that the saddle point 
corresponding to finite distance D approched the bulk saddle point when D — > oo. Our final 
results, Eqs ( 15T1 1521) confirmed the appearance of power law decay (~ D~ 3 ) of the Casimir 
force in the presence of condensate (/i > an c ), and of exponential decay in the one phase 
region (/i < an c ). Here our predictions reproduced exactly the situation in a perfect Bose 
gas. However, we noticed a difference in the value of exponent v governing the divergence 
of the characteristic length scale at the approach from below to the critical value of the 
chemical potential, /i c = an c . In a perfect gas vp = 1/2, whereas in the presence of the 
repulsive mean- field vimp = 1, as shown in Eq. (1481) . This difference seems to be the only 
manifestation of the mean-field coupling between the bosons as far as the Casimir effect is 
concerned. 

Finally, let us make the following comment. Because of non-equivalence of the grand 
canonical and canonical ensembles of a perfect gas one could wonder how the results derived 
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in Ref.- would compare with calculations based on the canonical free energy. It turns out 
that the answer to this important question can be readily derived from the remarkable, 
rigorous analysis of the properties of a perfect Bose gas in a finite volume given in Ref.—, 
see especially Section 3. We checked that as far as the Casimir force is concerned both grand 
canonical and canonical ensemble give exactly the same results. The anomalies of the grand 
canonical fluctuations in the total occupation number of the ground state in the presence of 
condensate have thus no effect on the large distance behavior of the Casimir force. 
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